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Abstract. Mass action systems capture chemical reaction networks in homogeneous and dilute 
solutions. We suggest a notion of generalized mass action systems that admits arbitrary power-law 
rate functions and serves as a more realistic model for reaction networks in intracellular environments. 
In addition to the complexes of a network and the related stoichiometric subspace, we introduce 
corresponding kinetic complexes, which represent the exponents in the rate functions and determine 
the kinetic-order subspace. We show that several results of Chemical Reaction Network Theory carry 
over to the case of generalized mass action kinetics. Our main result essentially states that, if the 
sign vectors of the stoichiometric and kinetic-order subspace coincide, there exists a unique complex 
balancing equilibrium in every stoichiometric compatibility class. However, in contrast to classical 
mass action systems, multiple complex balancing equilibria in one stoichiometric compatibility class 
are possible in general. 

Key words, chemical reaction network theory, generalized mass action kinetics, complex bal- 
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1. Introduction. Dynamical systems arising from chemical reaction networks 
with mass action kinetics are the subject of Chemical Reaction Network Theory 
(CRNT), which was initiated by the work of Horn, Jackson, and Feinberg, cf. [2511211 
113] . In particular, this theory provides results about existence, uniqueness, and sta- 
bility of equilibria independently of rate constants (and initial conditions). However, 
the validity of the underlying mass action law is limited; it only holds for elementary 
reactions in homogeneous and dilute solutions. In intracellular environments, which 
are highly structured and characterized by macromolecular crowding, the rate law has 
to be modified, cf. 0[2j[28]. 

Two types of modifications have been proposed: "fractal reaction kinetics" [251 
HHEElin] and the "power-law formalism" [MJ [35J |M1 EH ■ The names of the two 
approaches are a bit misleading, since both approaches address the problem of di- 
mensional restriction (i.e. molecules confined to surfaces, channels, or fractal-like 
structures) and both use power-laws. More specifically, in fractal-like kinetics, rate 
constants are time-dependent (via a power-law) , whereas the exponents of the species 
concentrations in the rate function are the corresponding stoichiometric coefficients 
(as in mass action kinetics). On the other hand, in the power- law formalism, rate 
constants are time- independent (as in mass action kinetics), whereas the exponents 
of the species concentrations may be (nonncgative) real numbers different from the 
respective stoichiometric coefficients. For model selection, data have to be collected 
for many molecules and intracellular environments. Recent data of binding kinetics 
in crowded media [SJ [30] suggest that the power-law formalism is the preferred model. 

In this work, we study the consequences of the power-law formalism for chemical 
reaction networks. In particular, we demonstrate that several fundamental results 
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of CRNT carry over to the case of generalized mass action kinetics (i.e. power-law 
rate functions). There has been an early approach to account for generalized mass 
action kinetics [35], which entails a redefinition of the complexes of a network. Here, 
we suggest a different approach, where we keep the original complexes, but introduce 
additional "kinetic complexes" , which determine the exponents of the species concen- 
trations in the rate functions. This has the advantage that the underlying chemical 
reaction network and thus properties like weak reversibility and deficiency remain the 
same. 

From the kinetic complexes, we obtain (in addition to the stoichiometric subspace) 
a "kinetic-order subspace" , and it turns out that the generalization of a central re- 
sult of CRNT (concerned with the uniqueness and existence of a complex balancing 
equilibrium in every stoichiometric compatibility class) depends on the sign vectors 
of the two subspaces. Our main result Theorem 13.101 essentially states that, if these 
sign vectors are equal, there exists a unique complex balancing equilibrium in every 
stoichiometric compatibility class. In general, however, there may be more than one 
complex balancing equilibrium in a stoichiometric compatibility class, see Proposi- 
tion [XI] and Example 14.21 

Chemical reaction networks with non mass action kinetics are also studied in 
[SJ2JGI]. In this approach, one is interested in conditions that guarantee the uniqueness 
of equilibria. If autocatalytic reactions are excluded and if the dependence of the 
rate functions on the species concentrations corresponds to the stoichiometric matrix, 
the structure of the stoichiometric matrix alone guarantees uniqueness. Moreover, 
the properties of the stoichiometric matrix can be translated into conditions for the 
species reaction graph. As a consequence, this theory is applicable to many types 
of kinetics, however, it does not address the existence of equilibria. Existence and 
uniqueness of equilibria for general kinetics are discussed in [T2]. The methods are 
based on homotopy invariance of the Brouwcr degree in a related way to the approach 
in Section l3~3l 

Organization of the work. In the next section, we recall the definition of 
mass action systems and several fundamental results of CRNT. Then we introduce 
generalized mass action systems and discuss the results that carry over easily to this 
framework. In Section [31 we study uniqueness and existence of complex balancing 
equilibria; more specifically, we reformulate the problem and study injectivity and 
surjectivity of a certain map, (a simplified version of) which appears for example in 
toric and computational geometry or statistics. In Section|4j we discuss two examples 
of generalized mass action systems. Finally, we draw our conclusions and give an 
outlook to further lines of research. In the Appendix, we recall the relevant results 
on sign vectors of vector spaces and face lattices of polyhedral cones and polytopes. 

Notation. We denote the positive real numbers by M> and the nonncgative real 
numbers by M> . For a finite index set /, we write for the real vector space of formal 
sums x = Y^iEi x i ' w hh x, £ K, and R> and R> for the corresponding subsets. Given 
x e R 1 , we write x > if x <E K> and x > if x G M>. Further, we define e x e R> 
and hx(x) <E M. 1 componentwise, i.e. (e x )i = e Xi and (ln(x)), — hi(xi), the latter for 
x 6 Mi,. Finally, we define x o y £ M. 1 for x, y £ M. 1 as (x o y)j = x^ and x y £ R> for 
x, y £ R> as x y = Yl ieI x Vi , where we set 0° = 1. 
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2. Chemical reaction networks. In our presentation of CRNT, we follow the 
surveys by Feinberg (Ml [15] [16] and Gunawardena [22] . 

Definition 2.1. A chemical reaction network (^,^£,8^) consists of three finite 
sets: (i) a set 5? of species, (ii) a set^ C R>" of complexes, and (Hi) a set 8? <Z c €y. c € 
of reactions with the following properties: (a) Vy G : 3t/ G ^ smc/i i/iai {y,y') G ^ 
or (y', y) G ^ and f&j Vy e : (y, y) £ ^. 

Complexes are formal sums of species; they are the left-hand sides and right-hand 
sides of chemical reactions. For y G c & , we may write y = J2 s e.y V s s ' wnere Us is the 
stoichiometric coefficient of species s. As usual in chemistry, we write y —> y' for a 
reaction (y, y') G 8%. In a chemical reaction network, each complex appears in at least 
one reaction; moreover, there are no reactions of the form y —> y. 

A chemical reaction network ^ , 8$) gives rise to a directed graph with com- 
plexes as nodes and reactions as edges. Connected components L\, . . . ,L\ C ^ are 
called linkage classes, strongly connected components arc called strong linkage classes, 
and strongly connected components without outgoing edges T\ , . . . , T t C ^ are called 
terminal strong linkage classes. Each linkage class must contain at least one terminal 
strong linkage classes, i.e. t > I. The network {5^,^,8?) is called weakly reversible, if 
the linkage classes coincide with the strong linkage classes and hence with the terminal 
strong linkage classes. 

From a dynamic point of view, each reaction y —> y' G 8? causes a change in species 
concentrations proportional to y' — y G R" 5 ". The change caused by all reactions lies 
in a subspace of R 1 ^ such that any trajectory in R>" lies in a coset of this subspace. 

Definition 2.2. Let {5? \8%) be a chemical reaction network. The stoichio- 
metric subspace is defined as 

S = span{y' - y G R y ' | y ->■ y' G 81} . 

Further, let c' G R^f. The corresponding stoichiometric compatibility class is defined 
as 

(c' + S)> = (c' + S)nRf '. 



2.1. Mass action systems. The rate of a reaction y — > y' G ffl depends on 
the concentrations of the species involved. The explicit form of the rate function 
Jfcy^,y> '■ RjT — > R> is determined by the underlying kinetics. In the case of mass 
action kinetics, it is a monomial in the concentrations c G M>" of reactant species, i.e. 
J£y^ y i(c) = ky^yi c v with rate constant k y ^ y i G R>. In other words, the stoichiomet- 
ric coefficient of a species on the left-hand side of the reaction equals the exponent of 
the corresponding concentration in the rate function. It remains to formally introduce 
the rate constants. 

Definition 2.3. A mass action system (.y,'To,8l,k) is a chemical reaction net- 
work (5? ' ,^ 1 ,8%) together with a vector k G R'Sf of rate constants. 

Definition 2.4. The ordinary differential equation (ODE) associated with a 
mass action system (J^, 8%, k) is defined as 



with the species formation rate 

y^y'am 
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In order to rewrite the species formation rate, we use the unit vectors w y £ 
corresponding to complexes y £ c to and define 

• a linear map3 Y : M.^ -» M.^ with Yu y = y, 

• a nonlinear map 

• a linear map0 A: — > , x \-> ^y^y' x y (^V — 
Now, the species formation rate can be decomposed as 

r(c)= ky^y'^iv 1 -V) (2-1) 

= Y Yl k y^y' ° v ^y' ~ u y^ 

y->y'£!% 

= YA^(c) . 

Equilibria of the ODE associated with a mass action system satisfying A^f(c) = 
and c > are called complex balancing equilibria. The possibility of other (positive) 
equilibria suggests the definition of the deficiency of a mass action system. 

Definition 2.5. Let (J? ^ ,M,k) be a mass action system. The set of complex 
balancing equilibria is defined as 

Z={c£Rf \A$(c) = 0}. 

The deficiency of the system is defined as 

5 = dim(kcr(y) n im(^4)) . 



Originally, the deficiency was defined differently. As we will see in Proposition 
12.81 the two definitions coincide under certain conditions on the network structure. 
In Fig. 12.11 we summarize the definitions associated with a mass action system and 
depict their dependencies. 

Results. Now we are in a position to present several results of CRNT related to 
the Deficiency Zero Theorem. (The results are due to Horn, Jackson, and Feinberg 
[5S1[5U[T3]. For proofs, we refer the reader to the surveys [T3J[TII] or [12] •) As we will 
see later, corresponding statements also hold in the case of generalized mass action 
kinetics. We start with a foundational linear algebra result, which can be proved 
using the Perron-Frobcnius Theorem. 



1 The corresponding matrix amounts to Y sy = y s . 

2 The corresponding matrix amounts to A yy i = K y / y — S yy i 5Z !/ "g t if ^yy"' where K £ R* x 
with K yy i = ky^yi if y — » y' G ^ and K yy i = otherwise. 
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Fig. 2.1. T/ie mass action system {^,^,31, k): Associated 
(Definitions at arrowheads depend on tails.) 



and their dependencies. 



Theorem 2.6. Let (5^ , ^6 ' ,8%, k) be a mass action system with the associated map 
A, and let T\ , . . . , Tt C be the terminal strong linkage classes. Then: 

1. for i = 1, . . . ,t: 3xi & K> with supp(x^) = TJ 

2. kcr(A) =span{xi,...,xJ 

3. dim(ker(^)) = t 

The next result is an immediate consequence of Theorem 12. 61 

Corollary 2.7. Let (S fi , c £, 8$) be a chemical reaction network. If there exist rate 

constants k such that the mass action system (S^ , ^ , 2$, k) has a complex balancing 

equilibrium, then (S^,^^) is weakly reversible. 

If each linkage class contains exactly one terminal strong linkage class, the de- 
ficiency is independent of the rate constants and can be computed from basic pa- 
rameters of the chemical reaction network. The resulting formula was the original 
definition of the deficiency. 

Proposition 2.8. If a chemical reaction network (J? 9 , , 8&) is weakly reversible 
( or more generally if t = I), then, for all rate constants k, the deficiency of the mass 
action system {,5^ , %f \8$, k) is given by 5 = m — I — s, where m is the number of com- 
plexes, I is the number of linkage classes, and s is the dimension of the stoichiometric 
subspace. 

In the case of deficiency zero, weak reversibility guarantees the existence of com- 
plex balancing equilibria. 

Proposition 2.9. If a chemical reaction network (y is weakly reversible 

and S = 0, then, for all rate constants k, the mass action system (j^, '0> , ffl, k) has a 
complex balancing equilibrium. 

Theorem 12.61 further implies that the set of complex balancing equilibria can be 
parametrized by the orthogonal of the stoichiometric subspace. 
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Proposition 2.10. Let (^','Tf,M, k) be a mass action system with nonempty set 
Z of complex balancing equilibria. Then 

Z = {ce Mf | ln(c) - ln(c*) e S^} = {c*oe v \ ve S^} 

for any c* G Z . 

Finally, we recall a result concerned with the existence and uniqueness of a com- 
plex balancing equilibrium in every stoichiometric compatibility class. It can be 
proved using methods from convex analysis. 

Theorem 2.11. Let (,y,'^',M,k) be a mass action system with nonempty set 
Z of complex balancing equilibria. Then Z meets every stoichiometric compatibility 
class in exactly one point. 

In Section [3j we study the conditions under which a result analogous to Theorem 
12. Ill holds in the case of generalized mass action kinetics. 

2.2. Generalized mass action systems. Chemical reactions occur between 
entire molecules such that the stoichiometric coefficients are integers. Under the as- 
sumption of mass action kinetics, the rate functions are monomials in the concentra- 
tions of the reactant species. However, in Definition 12 . 1 1 we allowed nonncgative real 
stoichiometric coefficients and hence "generalized monomials" as rate functions, since 
all results presented above also hold in this generality. This observation can be used 
to account for generalized mass action kinetics. We outline two different approaches 
the second of which is the focus of this paper. 

In the first approach |25) . chemical reactions are redefined as pseudo- reactions 
with the same net balance, but real stoichiometric coefficients. For example, the 
reaction 

uaA + nsB — > ncC 

with riA, riB, n>c 6N can be redefined as 

v A A + v B B + v c C ->■ (v A - n A )A + (v B - n B )B + (v c + n c )C 

with va,v Bi vc E K> and rate function k [A] UA [B] VB [C] uc . The redefinition of chem- 
ical reactions does not affect the stoichiometric subspace, however, it entails a new 
(and typically larger) set of complexes and hence a new mass action system (with 
different properties). For example, consider the (weakly) reversible chemical reaction 
network 

A + B^C 

with two complexes and one linkage class. Since the stoichiometric subspace S = 
span{(— 1, —1, 1) T } has dimension one, we obtain 5 = 2 — 1 — 1 = by Proposition ^. 81 
In order to account for generalized mass action kinetics specified by the rate functions 
kA+B^c[A] a [B] b and fcc^A+s[C] c with a, b 7 c € K>, the system can be redefined by 
the pseudo-reactions 

aA + bB (a - I) A + (b - 1)5 + C 
cC -+A + B+(c-l)C 

with four complexes and two linkage classes. This new system is not weakly reversible 
and has deficiency (5 = 4 — 2 — 1 = 1, again by Proposition 12.81 
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In this paper, we present a different way to account for generalized mass action 
kinetics. Most importantly, we disentangle the definition of the rate functions from 
the stoichiometric coefficients. In particular, we keep the integer stoichiometric coeffi- 
cients, but we allow "generalized monomials" as rate functions, in which the exponents 
of the concentrations can be arbitrary nonncgative real numbers. More formally, we 
do not change the chemical reaction network, but we associate with each complex a 
so-called kinetic complex, which determines the exponents of the concentrations in 
the rate function of the respective reaction. In the above example, we associate the 
kinetic complexes a A + bB and cC with A + B and C, thereby specifying the rate 
functions kA+B^c[A] a [B] b and fcc-^A+s[C] c . We obtain the following network, where 
we indicate association of kinetic complexes by dots: 

A + B ^ C 
aA + bB cC 

For an arbitrary chemical reaction network with generalized mass action kinetics, the 
rate function Jffy^y' ■ R>" — > R> corresponding to reaction y — > y' £ 3i is given by 
J£y^ y '(c) = ky^yi c v , where y is the kinetic complex associated with y. 

Definition 2.12. A generalized chemical reaction network (5^ , c tg, 38) is a 
chemical reaction network {^,^,38) together with a family c (o = (x v ) y ^ in 

Rf of 

kinetic complexes, where \{x y \y £ ^}| = \ c (?\. We write y = x y for the kinetic 
complex associated with the complex y £ c & . 

A generalized chemical reaction network {5? ', ^ , ^ 38) contains the chemical reac- 
tion network (5^,^,38); moreover, it entails the fictitious chemical reaction network 
(5?, , 3£) where the set ^? = {y | y £ has the same cardinality as ^ (by defini- 
tion) and the relation 3?_ is isomorphic to ffl, i.e. y — > y 1 £ M_ whenever y — > y' £ 3%. 
Hence the networks (=5^,^,^) and (=5^,^,^) give rise to the same directed graph 
(up to renaming of vertices). A generalized chemical reaction network (5^ ^ ^ , 38) 
is called weakly reversible if (^, ^ ', 38) is weakly reversible. Also the definitions of the 
stoichiometric subspace and the stoichiometric compatibility classes carry over from 
{^,^,38) to {^,^,^,38), cf. Definition O Additionally, we introduce the kinetic- 
order subspace of a generalized chemical reaction network, which coincides with the 
stoichiometric subspace of the fictitious network. 

Definition 2.13. Let (S* , < €, 3%) be a generalized chemical reaction network. 
The kinetic-order subspace is defined as 

S = span{y — y\y — s> y' £ 38} . 

For consistency, the name kinetic subspace would be more appropriate for S but 
this name has already been given to a certain subspace of the stoichiometric subspace 
[TB] . which coincides with the stoichiometric subspace if t = I. 

For later use, we introduce the maps 

• Y : M.^ -> with Yw y = y and 

• * : chj^c 4 uj y , 

y€l? 

where we identify R"^ and E— . 

Definition 2.14. A generalized mass action system {5? , 3%, k) is a gener- 
alized chemical reaction network (S^, 3i) together with a vector k £ R^ of rate 
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constants. 

Definition 2.15. The ordinary differential equation (ODE) associated with a 
generalized mass action system (S^, , k) is defined as 



dc 

dt = r(c) 



with the species formation rate 



As in Eqn. (|2.ip . we can decompose the species formation rate of a generalized 
mass action system as 

f(c) = YA^(c) . 

In analogy to Definition 12. 51 equilibria satisfying A$(c) = and c > are called 
complex balancing equilibria; they coincide with the complex balancing equilibria of 
the fictitious mass action system 3%_, k). The deficiency, which quantifies the 

possibility of other equilibria, coincides with the deficiency of the mass action system 
(J^,%f,^,fc). 

Definition 2.16. Let (.y,^,^,^!^) be a generalized mass action system. The 
set of complex balancing equilibria is defined as 

Z = {cGKf \A$(c) = 0} 

and the deficiency as 

5 = dim(ker(F) n im(A)) . 



It remains to introduce the kinetic deficiency, which coincides with the deficiency 
of the fictitious system. 

Definition 2.17. Let , 3$ , k) be a generalized mass action system. The 

kinetic deficiency is defined as 

8 = dim(ker(Y) n im(A)) . 



In Fig. 12. 21 we summarize the definitions associated with a generalized mass action 
system and depict their dependencies. From the mass action system (S^, 'if , 30, k), we 
keep the stoichiometric subspace S and the deficiency 8, whereas we use all defini- 
tions associated with the fictitious mass action system &)i m particular, the 
kinetic-order subspace S, the kinetic deficiency 8, and the set Z of complex balancing 
equilibria. 

Results. Now we return to the results of CRNT that have been derived for mass 
action systems. Since Theorem 12.61 is concerned with the kernel of the linear map 
A, the underlying kinetics is not relevant at all. But also Corollary 12.71 and Proposi- 
tions I2.8H2.101 carry over easily to generalized mass action systems if we consider the 
fictitious chemical reaction network M) and the fictitious mass action system 

k) defined above. For reference, we present the analogous results. 
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Proposition 2.18. Let (y,^,S?) be a chemical reaction network. If there exists 
a generalized mass action system (y, Si, k) with a complex balancing equilib- 
rium, then {y , ^ , Si) is weakly reversible. 

Proof. Assume that (y, ^ ', 'tf, Si, k) and hence the mass action system {y , Si_, k) 
have a complex balancing equilibrium. By Corollary 12.71 the chemical reaction net- 
work [y , c £_. Si) and hence (y, c £ , Si) are weakly reversible. □ 

Proposition 2.19. If a chemical reaction network {y ^ ,Si) is weakly reversible 
( or more generally ift = l), then the deficiencies of any generalized mass action system 
(y, ^ ', ^ , Si, k) are given by 5 = m — l — s and 5 = m — l — s, where m is the number of 
complexes, I is the number of linkage classes, s is the dimension of the stoichiometric 
subspace, and s is the dimension of the kinetic-order subspace. 

Proof. Assume that (y, ^ , Si) and hence the chemical reaction network (y , M ) 
arising from (y, ^ , ^ ,Si, k) are weakly reversible (or more generally that t = I). The 
deficiency of the generalized mass action system equals the deficiency of (y, 'Hi ' ,Si,k), 
and the kinetic deficiency equals the deficiency of (y, k). By Proposition 12.81 

the deficiencies of the two mass action systems are given by the formulas stated. □ 

PROPOSITION 2.20. If a generalized chemical reaction network (y, c ^,Si) is 
weakly reversible and 5 = 0, then any generalized mass action system (y, ^ , % ', SI, k) 
has a complex balancing equilibrium. 

Proof. Assume that (y, % Si) and hence the chemical reaction network (y,^_, S$ ) 
arising from (y, ff, ff, Si, k) are weakly reversible. Additionally assume (5 = 0. By 
Proposition [231 the mass action system (y, c £_, Si_, k) and hence (y, , ^ , Si, k) have 
a complex balancing equilibrium. □ 



10 



STEFAN MULLER AND GEORG REGENSBURGER 



Proposition 2.21. Let k) be a generalized mass action system with 

nonempty set Z of complex balancing equilibria. Then 

Z = {ce Kf | ln(c) - ln(c*) G S^} = {c*oe i |ie S^} 

for any c* G Z . 

Proof. The complex balancing equilibria of the mass action system ffl^ k) 

coincides with Z, and its stoichiometric subspace coincides with S, the kinetic-order 
subspace of 'to ', "W, k). By Proposition ^. 101 the nonempty set Z is given by the 
formula stated. □ 

One might conjecture that also Theorem 12.111 holds for generalized mass action 
systems. However, an analogous result depends on both the complexes and the 
kinetic complexes where C H> determines the stoichiometric subspace S (and hence 
the stoichiometric compatibility classes (d + S)>), whereas determines the set Z 
of complex balancing equilibria (and the related kinetic-order subspace S). It turns 
out that the result depends on additional assumptions concerning the sign vectors of 
the subspaces S and S, see Theorem 13. 101 

3. Complex balancing equilibria. In the following, we consider a generalized 
mass action system (S?, if, ', k) with stoichiometric subspace S, kinetic-order sub- 
space S, and nonempty set Z of complex balancing equilibria. 

From Proposition 12.211 we know that Z — {c* o e" | v G S^} for any c* G Z. We 
provide necessary and sufficient conditions such that in every stoichiometric compat- 
ibility class (d + S)> there is at most one complex balancing equilibrium. Moreover, 
we provide sufficient conditions such that in every stoichiometric compatibility class 
there is at least one complex balancing equilibrium. 

The question of uniqueness is answered by the following result for arbitrary sub- 
spaces S and S. It involves the corresponding sets of sign vectors denoted by <j(S) and 
o~(S); for the definition of sign vectors and related notions we refer the reader to the 
Appendix. We note that sign vectors also appear in the study of multiple equilibria 
that are not necessarily complex balancing (TTJ • 

PROPOSITION 3.1. Let S,S be subspaces o/R n . Then the two statements are 
equivalent: 

1. For all c* > and d > 0, the intersection (d + S)> fl {c* oe° \v G S- 1 } 
contains at most one element. 

2. <t(S) RfffS 1 ) = {0}. 

Proof, (-il => -i2): Suppose there exist b 1 / a 2 £ S and V\ ^ r>2 G S such that 
d + u 1 = c* o e" 1 and d + u 2 = c* o e V2 (for a certain d and a certain c*). Then 
u 1 — u 2 = c* o (e Vl — e" 2 ) and by the monotonicity of the exponential function 

aiu 1 - u 2 ) = a(c* o (e 5 ' - e" 2 )) = ^(e* 1 - e 5 *) = a(& - i 2 ) . 

Hence a(S) n a{S^) ^ {{)}. 

(-.2 => -.1): Suppose that / t £ ct(5) n c^S- 1 ). Then there exist it G S and 
w 1 G 5 X such that a(u) = ^(v 1 ) = r. Further, let v 2 = iv 1 . Then a(v x —v 2 ) = t and 

a{u) = aiv 1 - i 2 ) = a(e cl - e* 2 ) = a(c* o (e 5 ' - e fi2 )) 

for all c* > 0. In particular, there is c* such that u = c* o (e yl — e" 2 ). With d = 
c* o e" , one has d — u = c* o e" and hence both d and d — u arc elements of 
(c' + 5)>n{c*oc' D I^gS 1 - 1 }. □ 
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It follows in particular that if the sign vectors are equal, cr(S) = cr(S), complex 
balancing equilibria are unique (in a stoichiometric compatibility class) since then 

<t(S) n a(S^) = <t(S) n a(S) 1 - = <t{S) n a(S) 1 - = {0} 

using Eqn. (|A.1[) . Note that this is only a sufficient condition; for example, with 
S = span{(-l, 1)} and S = span{(-l, 0)}, we have a(S) n ^(S*- 1 ) = {0} but a(S) ^ 
cr(S). However, it includes classical mass action kinetics where S = S and each 
stoichiometric compatibility class contains at most one complex balancing equilibrium. 
On the other hand, if cr(S) D o-(S J -) ^ {0} and the underlying network is weakly 
reversible, then such a generalized chemical reaction network has the capacity for 
multiple complex balancing equilibria, as shown in the following result. 

Proposition 3.2. If a generalized chemical reaction network (J?, ff, c to, 3$) is 
weakly reversible and cr(S) D cr(S ) ^ {0}, there exist rate constants k such that the 
generalized mass action system (y, e ^ , , e if,I%,k) has more than one complex balancing 
equilibrium in some stoichiometric compatibility class. 

Proof. Let a(S) D ct^) ^ {0}. By Proposition O there exist c* > and 
d > such that (c' + S)> D {c* oe" | v S S^} contains more than one element. Using 
ProDOsition l2.21l it remains to show that there exist rate constants k <E R> such that 
c* is a complex balancing equilibrium of (J/*, "jf, ff, k), i.e. 

A*(c*) = ]T k y ^ y ,{c*)y{w y ,-w y ) = 0. 

Since {S?,^, and hence (J? ' ',&) are weakly reversible, this is guaranteed by 

Lemma l3~3l □ 

In the proof of Proposition ^. 21 we use the following result. 

Lemma 3.3. Let [5? ,M) be a chemical reaction network. Then, the following 
statements are equivalent: 

1. (.y ,^,1%) is weakly reversible. 

2. There exists k £ R> such that Y^y^y'zm ky->y' ( u y' ~ u y) = 0j where u y € R^ 
denotes the unit vector corresponding to y £ c £ ' . 

Proof. (1 => 2): By weak reversibility, there exists a cycle y — > y' — > . . . — > y for 
each reaction y —> y' 6 3% and we denote the set of reactions involved in this cycle by 
Cy^yi. Clearly, Sz->z'eC _> , ( w z' - w *) = and hence 

where k y ^ y i > records in how many cycles the reaction y — > y' appears. 

(2 => 1): We write Y, y ^ y >em k y^y' i^y'^^y) = A ^ with ^ = (1, 1, . . . , 1) T € K>. 
By Theorem if Aft = 0, then is weakly reversible. □ 

The second implication is a basic fact from CRNT [TB] . 

3.1. The map F. In order to study uniqueness and existence in a common 
framework, we rephrase the problem. We suppose that 5? contains n species and fix an 
order among them. Then we can identify R r with R" such that S, S C R™. Further, 
let V = (v 1 , . . . ,v d ) and V = (v 1 , . . . , v d ) be bases for S 1 - and S ± , respectively. 
In other words, S 1 - = im(V) and dim(»S ) = d and analogously S = im(t^) and 
dim(5 ± )=<i. 
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An element in (d + S)>C\ {c* a e v \ v £ S- 1 } corresponds to u £ S and v € 5 
such that c* o e" = c' + u or equivalently to A G R d such that 

(c'oc^^^Xc',!!') for i = l,...,d. 

Hence, provided c* e uniqueness and existence (of a complex balancing equilibrium 
in every stoichiometric compatibility class) correspond to injectivity and surjectivity 
of the following map: 

F: K d ~ ^ C° C K d (3.1) 
A h-» F(\) with (F(X)) t = (c* oc^U x ^\v l ) , 

where c* > and 

C° = {< 7 eR d \ li = (c',v i ), c'eR^}. 

Note that F depends on c*. It is instructive to reformulate the definition of F. To 
this end, we express the columns of V and V by its rows, 

V=(v\...,v d ) = (w\...,w n ) T 

v = (v\...y) = (w\...,w n ) T , 

or equivalently v\ — Wj and v\ = t/)*-, and obtain: 

(F{Xj) t = (c'oc^iV^'j = /Z c *k c^U x ^l v \ 

k=l 

n n 
fc=l A-=l 

and 

n n 

7* = (C , W ) = 2^ C fe W fc = 2^ C fe ™* ■ 
fc=l fc=l 

Hence we can write F{\) = J22=i c fc c^ A '"' 1 ^ w k and 7 = X)lLi c 'k wk ■ 

Definition 3.4. LetV EW lxd , V eR nxd with n> d,d have full rank. We write 
V = (v\...,v d ) = (w\...,w n ) T andV = (v 1 ,...,^) = (w 1 , . . . , w n ) T . Further, let 
c* > 0. We define 

F: R d -> C° C R d 

n 

where 

n 

C° = {^c' k w k G M d |c' 6 1",}. 
fc=i 
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This definition is more transparent than the equivalent one given above. It be- 
comes clear that the set C° is the interior of the polyhedral cone generated by the 
vectors (w 1 , ...,w n ). The map F itself (in case V = V) appears in toric geometry 
|20| . where it is related to moment maps, and in statistics [31] . where it is related to 
exponential families. There is a useful result [3D], which guarantees injectivity and 
surjectivity of F in case V = V. 

Proposition 3.5. Let V, V, and F be as in Definition \3.4\ IfV = V, then F 
is a real analytic isomorphism ofM. d onto C° for all c* > 0. 

This is a variant of Birch's Theorem [31] HUGS; it implies Theorem l2.11l We will 
build on this result when we study the surjectivity of F, but first we deal with its 
injectivity in case V ^ V. 

3.2. Injectivity of F. In the context of multiple equilibria in mass action sys- 
tems [TU| and geometric modeling [TT], it was shown that the map F (in case d = d) 
is injective for all c* if and only if F is a local isomorphism for all c*. We give an 
alternative proof of this result and extend it to the case d ^ d, where we use the sign 
vectors of the spaces im(F) and im(V). 

Theorem 3.6. Let V, V, and F be as in Definition \3.4\ Then, the following 
statements are equivalent: 

1. F is injective for all c* > 0. 

2. F is an immersion for all c* > 0. is injective for all A and c* > O.J 

3. cr(im(V) x ) n er(im(V)) = {0}. 

Proof. We use F in the form of Eqn. p.l[) . 

(1^3): By Proposition EHJ 
Using S 1 - = im(V) and S 1 - — im(V), the injectivity of F for all c* is equivalent to 
the existence of at most one element in (c' + S)> D {c* oc" \ v G S ± } for all c' and c* . 

(-i2 =>■ -i3): Suppose that ^ is not injective (for a certain c* and a certain A), 
i.e. there exists a nonzero A' G W 1 such that ^-X' = 0. Since 




v' 



this is equivalent to the existence of c > and v' € im(V) such that (cov\ v l ) — for 
i = 1, . . . , d, which in turn is equivalent to c o v' £ im(V)' L . Clearly cr(c o v') = o~(v') 
and hence c^im^)- 1 ) n a{im(V)) ^ {0}. 

(-.3 => -.2): Suppose that ^ r e c^im^)- 1 ) n a(im(V)). Then, there exist 
u G im(F) ± and v' G im(V) such that a(u) = cr(v') = r. Clearly, one can choose 
c > such that u = c o v' and hence c o v' G ixa(V) . As demonstrated in the 
previous step, this is equivalent to the existence of c* > and A, A' ^ G M d such 
that |f A' = 0. □ 

Finally, we note that for d = d, Statement 3 in Theorem 13.61 is symmetric with 
respect to V and V. 

Corollary 3.7. Let V , V be as in Definition \3.J\ with d = d. Then, a (im(V ") ± )n 
o-(im(V)) = {0} if and only i/a(im(y) ± ) n a(im(V)) = {0}. 

Proof. Let F be in the form of Eqn. (|3.ip with d — d, and let F be obtained from 
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F by changing the roles of V and V, 
F: R d R d 

X^F(X) with (F(X))i = (c* o e^i =1 Xjyi , v i ) . 

We will show that ^ is injective for all c* if and only if is injective for all c*. 
Then, by Theorem 13.61 we will obtain the desired result. 

Suppose that (or equivalently its transpose) is not injective (for a certain c* 
and a certain A), i.e. there exists A' £ K d such that (§j) T A' = 0. Since 

3=1 1 3=1 c 3=1 

and (c o t> 1 , {;'} = (c o v', v l ), this is equivalent to the non-injectivity condition for 
derived in the proof of Theorem 13.61 □ 

3.3. Surjectivity of F. It is more difficult to derive conditions for the surjec- 
tivity of F. Our main result is concerned with sufficient conditions, however, we start 
with a discussion of necessary conditions. 

Let C and C be the polyhedral cones generated by the vector configurations 
V T = (si) 1 ,...,™") and V T = (w 1 ,...,^), respectively. Then C° = int(C), and 
analogously we write C° = int(C). We note that C° and C° are nonempty since V 
and V have full rank. 

We write a(im(V))> = cr(im(y)) n {0, +} n for the face lattice of C, see the 
Appendix, and analogously cr(im(V"))> = a(im(V)) D {0, +}" for the face lattice of C. 
A face / of C is characterized by a sign vector r s <r(im(V^))> or equivalently by a 
supporting hyperplane with normal vector A £ M. d , where = whenever (A, w k ) = 
(for w k lying on /) and 7> = + whenever (A, w k ) > 0. 

Now, we can study a necessary condition for surjectivity: The image of F must 
contain points arbitrarily close to any point on a face of C. We assume that C is 
pointed, more specifically that (+, . . . , +) T € <j(im(V)), and we consider the simplest 
nontrivial face, namely an extreme ray e. To begin with, we assume that e contains 
only one generator, say w ; hence the characteristic sign vector amounts to r = 
(0, +, . . . , +) T . If .F is surjective, then the cone C must have a corresponding extreme 
ray e with the same sign vector r. Only then there is \x S K d with (/XjW 1 ) = and 
(/j,,w k ) > for = 2, . . . , n such that the limit 

n 

lim F(-ap + i/) = lim VcJ e - a ^ ,B ' = >+< ,/ ' ,s ' s > w fe = c\ e^'^ w 1 

a— >oo a— >oo ^ — ' 

fc=l 

can be placed arbitrarily close to any point on e (by appropriate choice of v £ 

If the extreme ray e contains more than one generator, there may be several 
corresponding extreme rays e. For a particular e with characteristic sign vector f, 
there is fj, £ R d (with (fi,w k ) = if ffc = and (/i,w fc ) > if ffe = +) such that 
lmio-^oo F{—a^L-\-v) lies on e. We note that if e contains then e must contain 
otherwise the limit does not lie on e. This condition on the corresponding extreme 
rays e and e can be expressed by their characteristic sign vectors f and r, namely as 
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f > r. For higher-dimensional faces, similar (but more complicated) conditions can 
be formulated. 

For the proof of the following surjcctivity result, we will employ Degree Theory. 
In particular, we use two properties of the Brouwer degree d(/, D, y) of a continuous 
function / : D — > M. d defined on the closure of an open and bounded subset D C R d 
(with boundary dD) at a value y g" f(dD): (i) the degree is invariant under homotopy, 
and (ii) if the degree is nonzero, there exists x such that y = f(x), see [29] or [19] . 

Theorem 3.8. Let V , V , and F be as in Definition \3-4\ If there exists a 
lattice isomorphism $ : cr(im(V))> — > a(im(V))> with f > $>(f) and (+,..., +) T G 
a(im(V)), then F is surjective for all c* > 0. 

Proof. In order to use the Brouwer degree, we require a map on a closed and 
bounded set. To this end, we define a map G equivalent to F from the interior of C 
to the interior of C and extend G to the boundaries such that it maps faces to faces. 
Then, we cut the pointed cones such that we obtain polytopes P and P. Finally, we 
define a homotopy between the map G and a homeomorphism between the polytopes 
guaranteed by the face lattice isomorphism. As a consequence, every point in the 
interior of P has nonzero Brouwer degree and hence is in the image of G. Since the 
cut of the cone C can be placed at arbitrary distance from the origin, this holds for 
every point in the interior of C . 

Since (+, . . . , +) T € er(im(F)), the face lattice isomorphism implies (+, . . . , +) T G 
a(im(V)), and hence the cones C and C are pointed. We start by choosing a minimal 
set of generators for C, which (after reordering) we assume to be (w 1 , . . . , w nE ), where 
he is the number of extreme rays of C. We define an auxiliary map, 

F: R d ~ -> C° 

n E 
k=l 

which is a real analytic isomorphism by Proposition 13.51 and a composed map, 

G°:C°^ C° 

F(F-\x)), 

which is surjective whenever F is surjective. 

Since G° is defined only on G° , we want to extend it continuously to the bound- 
ary dC, i.e. to the faces of the cone. Let / be a face of C. It contains a subset 
of the minimal set of generators for C, which (after reordering) we assume to be 
(w , . . . , ?D™ min ). There may be additional generators on /, which we assume to be 
(w nE+1 , . . . , u}"E+"add) ) where n min + rt a dd is the total number of generators on /. 

Now, let {x 1 )^ be a sequence with x % G C° and lim.^oo x % G /. Via the isomor- 
phism F, there is a corresponding sequence (A 1 )^ with A 1 G M. d . From 

lim x i = V c% lim c^^^ w k + V c* k lim e < A *'**> w k , 

fe=l /c=n min + l 

we conclude that lim^ooC^ ' w1 "> > for k — l,...,n m in and lim^oo e^ A ' w ) = 
for k = n m i n + 1,...,He- Additional generators w k on / can be written as non- 
negative linear combinations of the minimal generators (w 1 , . . . , w n ™ ia ) and henc^l 

3 By using c( A 'iT = ]J^f (e^\™ k )y k . 
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we obtain limj_}. 0o e^*'' 1 " ^ > 0. Generators w k not on / can be written as nonnega- 
tive linear combinations containing at least one of the remaining minimal generators 
(ui™ min+1 , . . . , w nE ) and henc^l we obtain limi_ ) . 00 e( y > wk ) =0. As a consequence, the 
image of the sequence converges and 

«min »E+"add 

lim G°(x i ) = V c* k lira e^'™^ w k + V c* k lim e {yx " k) w k . 

The isomorphism $ (between the face lattices of C and C) with f > ^(t) implies 
that there is a face / of C with w fe 6 / if u> fc £ /. That is, io 1 , . . . , iu™ min £ / as well 
as . . . , w ™B+™add e j and hence lim-^^ G° (x*) e /. 

In other words, there is a continuous extension of G° to the face /, which maps / 
to the corresponding face /. We set G :~ G° on C° and G(x) := lim^oo G°(x l ) for 
any sequence (x l )i^ with x % £ C° and lim^oo x l = x £ /. Since this can be done 
for all faces of C, there is a map G: C — » C which extends G° continuously to <9(7 
and maps faces to faces. 

Due to the face lattice isomorphism, a minimal set of generators for C is given 
by (w , . . . , w riE ). The isomorphism further implies d = d. 

Since C is a pointed cone, we can choose & (d — l)-dimcnsional subspace of M. d 
such that C lies on one side of the subspace. We cut C with a hyperplanc parallel 
to the subspace and obtain a polytope P (lying on one side of the hyperplanc). In 
particular, we intersect the extreme rays of C with the hypcrplane: the intersection 
of the extreme ray e k (generated by w k ) is located at atw k with au > 0. 

Analogously, we cut C with a hypcrplane and obtain a polytope P. The intersec- 
tion of the extreme ray e k (generated by w k ) is located at cekW k with au > 0. 

From now on, we restrict the map G to P and choose c* such that G maps 
corners of P to corresponding corners of P. For example, the corner ajiii 1 on e 1 
corresponds (by to the sequence (A l )i 6 N with A* £ R", lim^oo c* ' w ) = ai, 

and lim^oo e^ A ,w ^ = for fe = 2, . . . , Ue- In turn, (A l )i S N corresponds (by F) to the 
corner aiw 1 on e : 



it; | = ot\W 



lim cT e^'^ 1 ) + V ct c< v ^> «.* 1 = — 1 



Here, we have assumed that in addition to w 1 there are additional generators w k 
(with k = tie + 1, . . . , he + riadd) on e 1 with corresponding generators w k on e 1 . If 
we write w k = PkW 1 , w k = PkW 1 , and x = lim^oo e^ A ' w \ we can determine from 

n.E+™add 

x = 6t\ with x + c* k x^ k ftk = «i • 

k— riE + 1 

If we choose c£ accordingly for each extreme ray , then G maps "side-edges" of P 
to corresponding side-edges of P. The image of other faces of P need not coincide 
with the corresponding faces of P. (However, due to the face lattice isomorphism, 
the image of a "side-face" of P lying on a face of C, lies in the corresponding face of 



By using e<^Ki ^ = U^T (e<*'.**>)"* Yll^+x (e<*'.**>)°* 
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C.) In particular^, the image of the "cut-face" of P (arising from the cut with the 
hyperplane) may lie outside the cut-face of P. 

The isomorphism between the face lattices of C and C has another important con- 
sequence. It guarantees the existence of a piecewise linear homeomorphism G' : P — > 
P, which restricts to homeomorphisms between corresponding faces of P and P, see 
the Appendix. We note that G' has nonzero Brouwer degree on P° = int(P) and 
define a homotopy between G (restricted to P) and G' , 

H: P x [0,1] C C R d 

(x,t)^-tG(x) + (l-t) G'(x). 

(The homotopy H maps to C, since both G and G' map to C and C is convex.) 

Now, let y £ P°. Below we will show that y g H(dP,t) for all t G [0, 1]. Writing 
P° = int(P), we conclude that d(G,P°,y) = d(G',P°,y) ^ (by the homotopy 
invariance of the Brouwer degree) and that there exists x G P° with G{x) = y (by the 
existence property of the Brouwer degree). In other words, the image of G restricted 
to P° contains P°. Since the cut of the cone C can be placed at arbitrary distance 
from the origin, G° : C° — > C° and hence F : M. d — > C° are surjective. 

It remains to show that y ^ H(dP,t) for all t G [0,1]: For side-faces / C dP, 
one has H(f,t) C dC for all t 6 [0, 1] (since G and G' map side-faces to side- faces), 
whereas for the cut-face / C dP, one either has H(f,i) C dP for all t G [0,1] 
(whenever G maps one cut-face to the other) or H(mt(f), t) n P = for all t G [0, 1] 
(whenever G maps the cut- face out of P). In each case, one obtains H(dP, t)P\P° = 
for all t G [0,1]. □ 

We think that the technical condition (+, . . . , +) T G a(im(V)) in Theorem 13.81 
which requires the cone C to be pointed, is not necessary, and a similar result can 
be obtained for arbitrary cones. However, at the moment we do not have a complete 
proof for such a theorem. 

3.4. Main results. The previous two theorems concerned with injectivity and 
surjectivity of F allow the following generalization of Proposition 13.51 (Birch's Theo- 
rem). 

PROPOSITION 3.9. Let V, V, and F be as in Definition \3.4\ If a (im(V)) = 
<r(im(V')) and (+, . . . ,+) T G o~(im(V)), then F is a real analytic isomorphism ofM. d 
onto C° for all c* > 0. 

Proof. From <r(im(V)) = o-(im(V)) it follows that d = d and with Eqn. (|A.1[) 
that a(im(V) ± ) n cr(im(y)) = {0}. Hence, F is injective and a local isomorphism by 
Theorem 13.61 Moreover, with $ being the identity, F is surjective by Theorem l3.8l □ 

Note that the condition <r(ini(V)) = <r(im(V)) in the previous proposition can be 
tested algorithmically using chirotopes, see the Appendix. We can now formulate a 
result analogous to Theorem 12. Ill in the case of generalized mass action kinetics. 

Theorem 3.10. Let (5? , k) be a generalized mass action system with 

nonempty set Z of complex balancing equilibria, stoichiometric subspace S, and kinetic- 
order subspace S. If <r(S) = cr(S) and (+,..., +) T G a(S ), then Z meets every 
stoichiometric compatibility class in exactly one point. 



5 A point on the cut-face of P is a convex combination of the "corners" a^ui*, k = 1, . . . ,tie- 
By F~ 1 it corresponds to some A £ IR* 1 , which by F corresponds to a point on the cut-face of P, 

that is, a convex combination (with the same coefficients) of the corners a^w k , k = 1, , We, plus 

a positive linear combination of the additional generators w k , k = + 1, . . . , n. 
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Proof. Suppose Z ^ 0. As discussed at the beginning of Subsection 13.11 unique- 
ness and existence of a complex balancing equilibrium in every stoichiometric com- 
patibility class correspond to injectivity and surjectivity of the map F as given in 
Definition 13.41 where V and V are bases for S and S , respectively. By Eqn. (|A.1[) . 
a(S) = (t(S) is equivalent to a(im(V)) — er(im(V)), and obviously (+,..., +) T G 
o~(S ± ) is equivalent to (+, . . . , +) T £ cr(im(y)) such that F is injective and surjective 
by Proposition 13.91 □ 

In the terminology of CRNT, a chemical reaction network is conservative if S n 
^ 0, i.e. if there is a "vector of molecular weights" , relative to which all reactions 
are mass conserving. Note that the condition (+, . . . , +) T £ cr(5 f± ) in Theorem 13. 101 
means that the underlying chemical reaction network is conservative. 

4. Examples. We discuss two examples of generalized mass action systems. 
First, we continue the example of the generalized chemical reaction network intro- 
duced in Section [2721 



A + B =^ C 

• : (4- 1 ) 

aA + bB cC 

with a, b, c £ K>. The kinetic complexes aA + bB and cC (associated with the com- 
plexes A + B and C) determine the exponents in the rate functions fc J 4+B-s.c[^l] a [-B] 6 
and k c ^A+B[C] c . 

The network is (weakly) reversible and has 2 complexes and 1 linkage class. The 
stoichiometric and kinetic-order subspace amount to S — span{(— 1, — 1, 1) T } and 
S = span{(— ua, —vb,vc) T } with dimensions d = d = 1. By Proposition 12.191 S = 
5 = 2 — 1 — 1 = 0, and by Proposition ^. 201 Z ^ 0. Further, the sign vectors of S and 
S coincide, i.e. a(S) = a(S), and (1, 1, 2) T £ S ± , which implies (+, +, +) T £ a(S ± ). 
Hence, by Theorem l3.101 every stoichiometric compatibility class contains exactly one 
complex balancing equilibrium. 

In the rest of this section, we study an autocatalytic mechanism (for the overall 
reaction A + B =^ C) endowed with generalized mass action kinetics: 

A + 2B =± B + C 

: : (4.2) 

A + B 2B + C 

The kinetic complexes A + B and 2B + C (associated with the complexes A + 2B 
and B + C) determine the rate functions kA+2B^B+c[A][B] and kB+c^A+2B[B] 2 [C]. 
The particular kinetics may be unrealistic from a chemical point of view, however, 
it will serve to demonstrate how the conditions in Theorem 13.101 for existence and 
uniqueness of a complex balancing equilibrium (in every stoichiometric compatibility 
class) are violated. 

The network is weakly reversible, 6 = 3 = 0, and hence Z 0. In particular, 
the stoichiometric and kinetic-order subspace amount to S = span{(— 1, — 1, 1) T } and 
S = span{(— 1, 1, 1) T }. For the orthogonal complements S and S 1 - we choose the 
bases 



V = 
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The cones C and C generated by V = (w ,w ,w°) and V = (w < ,w ,w°) 



both coincide with R> : 



„c 



w B 



w A 







1 







1 



First, we address the question of existence. We observe that the cone C has an 
extreme ray generated by w A , whereas the cone C does not have a corresponding 
extreme ray generated by w A (since w A lies in the interior of C). As a consequence, 
the map F is not surjective for all c*, cf. the argument at the beginning of Subsection 
13.31 In other words, there may be a stoichiometric compatibility class that does not 
contain a complex balancing equilibrium. 

Now, we turn to the question of uniqueness. In order to employ Proposition 13.11 
or !3.2[ we determine cr(S) (~1 cr(5 J '). The sign vectors of S are (— , — , +) T , its inverse, 
and 0, whereas the sign vectors of <S' J ~ can be read off from the above figure: For every 
hyperplane of IR 2 , i.e. for every line through E R 2 , we check if w A , w B , and w c lie 
on the line or on its negative or positive side. We obtain 



a(S) = 



and a(S^) = 
0/ 



+ + - - 
- - - - 
+ + + + 



where we use matrix notation for sets of vectors and where we do not state vectors 
explicitly that are inverses of others. We find that cr(S) H cr(5 ± ) contains (— , — , +) T . 
Hence, by Proposition 13.21 there exist rate constants kA+2B->B+c an d fcs+c->/i+2B 
such that some stoichiometric compatibility class contains more than one complex 
balancing equilibrium. 

Due to the simplicity of the generalized mass action system, the equilibria of the 
associated ODE can be determined analytically. The equilibrium condition amounts 
to 

k A+2 B^B+c[A]{B] = k B +C^A+2B[B} 2 [C] , 

and since 5 = all equilibria are complex balancing. By using the conservation 
relations [A] + [C] = [A] + [C] = ^ac and [B] + [C] = [B] + [C] = ^bc and by 
writing K = kA+2B^B+c /kB+c^A+2B, we obtain a quadratic equation in [C], which 
can be solved as 



[C] = 




- K Y. AC 



Depending on the equilibrium constant K and the initial values T>ac and £_bc (which 
determine a stoichiometric compatibility class), the quadratic equation has 0, 1, or 2 
solutions with [C] > 0. If additionally [A] = Y, AC - [C] > and [B] = E B c - [C] > 0, 
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then ([A], [B], [C]) T is a complex balancing equilibrium. Obviously, a stoichiometric 
compatibility class contains 0, 1, or 2 complex balancing equilibria; it turns out that 
each case is realized. 

5. Conclusion. CRNT establishes intriguing results about the ODEs associated 
with mass action systems, in particular about the existence, uniqueness, and stability 
of equilibria. For application in molecular biology, however, one would like to have a 
framework that permits rate laws more general than mass action kinetics. 

In this paper we show that the suggested notion of generalized mass action sys- 
tems, which admits arbitrary nonncgative power-law rate functions, allows to gener- 
alize several results of CRNT. In particular, Theorem 13.101 essentially states that if 
the sign vectors of the stoichiometric and the kinetic-order subspace coincide, there 
exists a unique complex balancing equilibrium in every stoichiometric compatibility 
class. 

A natural next step is to study other results of CRNT in the case of generalized 
mass action kinetics, most importantly, to analyze the stability of complex balancing 
equilibria, which is guaranteed in the classical case. Further, genuinely biological no- 
tions such as the robustness [6l[39l|40] of chemical reaction networks can be addressed 
in a framework with more realistic kinetics. 

Appendix. Sign vectors and face lattices. 

In this section, we outline some facts on the relation between sign vectors of vector 
spaces and face lattices of polyhedral cones and polytopes. For further details we refer 
to fU Ch. 7] and |HJ Ch. 2, 6] and to [33J E] m the context of oriented matroids. 

We obtain the sign vector a(x) G {—,0, +}" of a vector x G R n by applying the 
sign function componentwise, and we write 

a(S) = {a(x) | x G S} 

for a subset SCI". 

Two sign vectors <r,r G { — ,0, +}" are orthogonal, if = for all k or if there 
exist k, I with Q-Tfe = — and qti = + (where the product on { — ,0,+} is defined in 
the obvious way); we write <r^r. Note that ^_Lr if and only if there arc orthogonal 
vectors x, y G W 1 such that a(x) = g and a(y) = r. 

The orthogonal complement E 1 - of a set E C {— , 0, +}" is defined by 

E- 1 = {<; G {-, 0, +}" | <r_Lr for all r G E} . 
The sign vectors of the orthogonal complement of a subspace S C R™ are given by 

o-^) = o-(S)^ (A.l) 

see for example [42j Prop. 6.8.]. 

Let V = (v 1 , . . . , v d ) G W ixd with n > d have full rank. Then V T = (u>\ . . . , w n ) 
is called a vector configuration (of n vectors in IR" 1 ). With A G R d and v = X)j=i ^j v ^ 
im(V), we obtain Vk = J2j=i ^j v l = Ej=i ^i w j = (^> wh )- Hence, a(v) describes the 
positions of the vectors w w n relative to the hyperplane with normal vector A. 

The face lattice of the cone C generated by w l , . . . , w n can be recovered from the 
sign vectors of the subspace generated by v 1 , . . . , v d . It is the set a(im(V)) fl {0, +}" 
with the partial order induced by the relation < +, which we denote by 

a(im(V))> = o-(im(V)) fl {0,+}". 
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A face / of C is characterized by a supporting hyperplane with normal vector A6R <i 
such that (A, w k ) = for generators w k lying on / and (A, w k ) > for the remaining 
w k (thus lying on the positive side of the hyperplane). 

A cone C is called pointed if C fl (— C) = {0} or cquivalcntly if it has vertex 0. A 
cone is pointed if and only if it has an extreme ray, and every pointed polyhedral cone is 
the conical hull of its finitely many extreme rays. Note that if (+, . . . , +) T <E <r(im(V)), 
the cone C generated by V T is pointed. 

As for polyhedral cones, the faces of a polytope form a lattice. Two polytopcs 
are combinatorially equivalent if their face lattices are isomorphic. Combinatorial 
equivalence corresponds to the existence of a piecewise linear homeomorphism between 
the polytopes that restricts to homeomorphisms between faces. 

The sign vectors er(im(V)) of the subspace im(V) can be equivalently charac- 
terized by the chirotope xv T °f the point configuration V T , which is defined as the 
map 

Xv t. {l,..., ? i} d ^{-,0,+} 

(ii, . . . ,i d ) i-> sign(det(w n , . . . , w ld )) . 

The chirotope records for each d-tuple of vectors if it forms a positively (or nega- 
tively) oriented basis of M. d or it is not a basis. It can for example be used to test 
algorithmically if the sign vectors of two subspaces are equal, that is, to decide if 
a(im(V)) = a{im(V)) for two matrices V, V € R nxd . 
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